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Abstract 

The unparticle is proposed by Georgi as a conceptual new physics beyond the 
standard model to describe the low energy sector of a nontrivial scale invariant 
sector of an effective theory. We consider the neutrino decay in unparticle scenario 
and investigate the effect of such decay for the case of ultra high energy neutrinos 
from GRB. The effect on the detection yield of these neutrinos are also probed for 
a kilometer scale ice-Cerenkov detector. 
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1 Introduction 



The scale invariance is a feature in which a theory or a law does not change 
if the energy scales or length scales are given a multiplicative transformation 
by a common factor. The transformation is also known as dilatation and 
this has applications in both mathematics and physics. In mathematics, 
the scale invariance can be described by self similarity. An object is such 
that if it is magnified by any amount, a smaller piece of the object always 
resembles the whole object. Fractals can be a typical example. In statistical 
mechanics, a scale invariant theory is needed to explain the phenomenon 
of phase transition as the fluctuations near the critical point become scale 
invariant. The scale invariance can apply in classical field theory, string 
theory etc. In quantum field theory, the scale invariance, can very naively 
be interpreted as the non-dependendence of the particle interaction strength 
on the energy of the particles. 

In the real world, the scale invariance is manifestly broken by the masses 
of the standard model particles. The scale transformation involves mass 
dimension and therefore in scale invariant scenario, the particle masses have 
to be zero. Therefore standard model does not have scale invariance. At 
a scale much higher than that of standard model, the scale invariance is 
restored. In a recent work [1, 2] Georgi observed that if an yet unseen scale 
invariant sector is present in four dimensions and very weakly interacting 
with the standard model particles, then this scale invariant sector contains 
not "particles" but massless "Unparticles" . 

Under the framework of field theory, an interacting scale invariant theory 
is invariant under conformal group transformation. A conformal group trans- 
formation is a simple transformation and can be generated (in Minkowsky 
space) by making a scale transformation of the type x^ 1 — > Aa; M (/x being the 
space-time coordinates) [3]. Although conformal invariance is well studied in 
two dimensional scenario, it is rare in four dimensional quantum field the- 
ory. The renormalization group effects break the conformal invariance [4]. 
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In field theory, renormalization group flows from some scale invariant ultra 
violet fixed points to some scale invariant infra red fixed points. However, 
this exceptional scenario of conformal invariance in four dimension can be 
described by the vector-like non abelian gauge theory with suitable number 
of massless fermions. Such a theory is first studied by Banks and Zaks (BZ) 
[5]. The theory has non-trivial infra red fixed points. At low energy limits 
a non-trivial conformal sector is ensured by the presence of non-trivial in- 
fra red fixed points in the theroy. But this requires the number of fermion 
generations to be non-integral number and therefore is not manifested in 
nature. 

As suggested in Ref. [1, 2] such a nontrivial scale invariant sector due to 
BZ field might appear at TeV scale. The BZ field interact with the standard 
model field via the exchange of a particle with large mass scale M u . This 
can be termed as connector sector [6]. Below M u standard model and BZ 
field couplings are suppressed by powers of M u . Generically the operators 
can be written as [1], OsmObz/M]), (k > 0) where Osm and Obz are op- 
erators of mass dimensions dsM and d^z respectively and constructed from 
standard model and BZ fields respectively. Below the connector sector such 
interactions between BZ and standard fields can be parametrized in terms 
of non-renormilazable interactions and the the suppression of nonrenormaliz- 
able operators by powers of Mjj are induced. As the scale invariance sets in at 
an energy scale A^, the renormalizable couplings of the BZ fields in the scale 
invariant sector induce dimensional transmutation [7] at A^. The particles in 
this sector is then described by massless unparticles [1]. In the effective the- 
ory, below the scale An, the BZ fields match onto the unparticle operators 
and a new set of operators having the form {Co u K dBZ ~ du / M^)OsmObz) , 
where du is the scale dimensions of unparticle operator On and Cq u is a 
coefficient to be fixed by the matching. This may be occured that the cou- 
pling between the unparticle and standard model particles may violate the 
conformal invariance but in [2] it is argued that the infrared fixed points of 
the unparticles will not be affected by such couplings since BZ fields decou- 



3 



pie from the standard model particles at low energies. As observed in Ref. 
[1], the unparticle stuff with scale domension du appears to be a nonintegral 
number du of massless particles. 

The unparticle physics opens up the possibility of very unexpected phe- 
nomenology and the effects of unparticle coupling for specific processes are 
studied by various authors [8]. Kikuchi and Okada [9] adressed the Higgs 
phenomenology with unparticles where they considered operators involving 
scalar unparticle, Higgs and the gauge bosons (gluons and photons) and dis- 
cussed the the effective couplings between the Higgs boson and the gauge 
bosons. The interaction of unparticles with Standard Model particles is also 
addressed by Chen and He [10]. Recently Zhou [11] discussed the possibiltiy 
of neutrino decaying into unparticle. Neutrino decay to unparticles are also 
discussed by Chen et al in Ref. [12]. If such neutrino decay really occurs it 
may have consequences in several ongoing neutrino experiments as also the 
future experiments. In a recent work Li et al [13] has discussed the neutrino 
decay in unparticle scenario as a possible explanation of the excess electron 
like events observed by the MiniBooNe [14] neutrino experiment. In the 
present work, the effect of such neutrino decay to unparticles is considered 
on ultra high energy (UHE) neutrinos arriving on earth from distant cosmic 
sources such as gamma ray bursts (GRBs). The decay length for decays 
to unparticles depends on quantities like the unparticle-neutrino coupling, 
the unparticle scale factor etc. The decay length for this scenario can be 
large enough (~ tens of Mpc) and thus can significantly affect the survival 
probability of a decaying neutrino if it traverses a baseline length of Mpc 
order. The ultra high energy neutrinos from distant GRBs at astronomical 
distances from the earth may indeed provide such a long baseline. Also, for 
such a long baseline the oscillatory effect of neutrinos are averaged out and 
therefore opens up the posibility to effectively probe such decays, if any, in 
signals from UHE neutrinos at terrestrial detectors such as IceCube [15]. 

The paper is organised as follows. In Section 2, the formalism is dis- 
cussed. This includes the neutrino flux from a GRB at a redshift z; the 
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survival probabilty of such neutrinos on reaching the earth if they undergo 
decay to unparticles. Therefore the actual flux of neutrinos on reaching the 
earth is obtained from folding the original flux with the survival probability. 
These are described in Section 2.1. In Section 2.2, the analytical expres- 
sions for the yield of secondary muons and shower events induced by the 
UHE neutrinos at the ice Cerenkov detector are described. The calculational 
procedure of the yield at the kilometer scale ice Cerenkov detector such as 
IceCube for the neutrino flux described in Section 2.1, is given in Section 
3. calculational results are also discussed in this section. Finally, Section 4 
contains discussions and summary. 

2 Formalism 

2.1 GRB neutrino flux with neutrino decay to unpar- 
ticles 

The GRB neutrino flux for a particular redshift is estimated considering the 
the relativistic fireball model [16]. In this type of GRB model, protons (also 
electrons, positrons and photons) produced in the magnetic field of the rotat- 
ing accretion disc around a possible black hole are accelerated perpendicular 
to the accretion disc at almost the speed of light. This forms a jet which is 
referred to as fireball. The burst is supposed to be the dissipation of kinetic 
energy of this relativistic expanding fireball. High energy pions are photopro- 
duced from A resonance when the protons in the jet interacts with photons. 
These pions then decay to yield and v e in the approximate proportion of 
2:1. 

The neutrino flux from a GRB depends on several GRB parameters like 
Lorentz boost factor T (required for the transformation from the fireball blob 
to observer's frame of reference), the photon break energy (as the photon 
spectrum is considered broken) and the photon luminosity L 1 (generally ~ 
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10 53 ergs/sec). 

With all these, the neutrino spectrum from a GRB can be parametrised 
as [17, 18] 

^ = Axmin(l,^)x-l. (1) 
In the above, E u is the neutrino energy N v is the number of neutrinos and 

El ~ 10 6 -p^GeV 

7,MeV 

r 2 . 5 = r/io 2 - 5 

= ^GRB , , 

1 + ln(^ max /^) ' 1 ; 

where E umSuX is the cut-off energy for the GRB neutrinos and -Egrb is the total 
energy that a GRB emits. Now, the observed energy E° bs of a neutrino with 
the actual energy E v coming from a GRB at a redshift distance z is given 
by the relation E° hs = E v /{1 + z) and similarly, the maximum observable 
neutrino energy E°^ x is E°^ x = E umax /(1 + z). The comoving distance d 
of a GRB at redshift z is given by 



<*(*) = ^ T / ^ 7 (3) 



where Om is the matter density of the universe at present epoch, Q\ is the 
dark energy density respectively in units of critical density of the universe 
and c,H are the velocity of light in vacuum and Hubble constant respectively. 
In the present calculation c = 3 x 10 5 Km/sec and H = 72 Km/sec/Mpc (1 
Mpc = 3.086 x 10 19 Km). Therefore the neutrinos from a single GRB that 
can be observed on earth per unit energy per unit area of the earth is given 
by, 

dN ohs dN 1 

ai\ v = aiv^ i , v 

dEf* dE v ^d\z) y } K } 

The production process of UHE neutrinos suggests that the neutrino 
flavours are produced in the ratio u e : : v T = 1 : 2 : 0. In case of 
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mass flavour oscillation, because of the astronomical baseline [L ~ Mpc), 
the acquired relative phases of the propagating neutrino mass eigenstates is 
averaged out (Am 2 L/E » 1) and the UHE neutrinos from a GRB reach- 
ing the earth are incoherent mixture of mass eigenstates. In fact, it can be 
shown [19] that for #23 = 45° (maximal mixing) and # 13 ~ 0, the flavour ratio 
on reaching the earth, for neutrino mass-flavour oscillation with such a long 
astronomical baseline, becomes v e : : v T — 1 : 1 : 1 irrespective of the solar 
mixing angle. 

However, in the present work we consider neutrino decay to the recently 
proposed "Unparticle" by Georgi and ivestigate that in the event of such 
decay, the possible signal of UHE neutrinos from distant GRBs in terms of 
muon and shower yields of a kilometer scale detector such as IceCube. 

We consider here a decay hypothesis where a neutrino mass eigen state 
Uj decays to an unparticle and another neutrino mass eigenstate z/j following 
the decay relation 

Vj + vi 

The unparticle physics proposed by Georgi [1] indicates a possibility for neu- 
trino decaying into the so called unparticle [11]. The effective Lagrangian for 
the process can be written as 

\<*P 

£ = -j^ZiValpVfiOu (5) 

where Ou is the unparticle operator, A„ is the relevant coupling constant and 
a, (3 are the flavour indices. In Eq. (7) dimension transmutation scale is An 
at which the scale invariance sets in and dy is the scaling dimension. 

In the mass basis, the interaction term (between neutrino and unparticle) 
can be written as 

\v 

jd^T^a^jOu (6) 
where i, j are the mass eigenstate indices and 

^ = E KX%> ■ (7) 

a,/3 
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In the above U, is the Maki-Nakagawa-Sakata (MNS) mixing matrix (mass- 
flavour) for neutrinos. The decay rate for the process Vj — > i/j + U (U is the 
unparticle) can now be elected [11] and the neutrino lifetime r for such decay 
process is given as 

ru_ = 2 d u ^d u {2-d u ){d u + l) ( A|\ du ~ X J_ 
3A d jA!/'| 2 \rnjj 

In Eq. (10) above, rrij is the mass of the decaying neutrino and A du is a 
normalization constant given by [1] 

lGvr 5 ^ rfrfo + 1/2) 
d " (2tt)^ r(rf w - l)V{2d u ) 1 J 

A decay scenario can be one in which both the states \u 2 ) and \u 3 ) are 
unstable and decay whereas the lightest state \ui) is stable. Needless to say 
that the normal mass hierarchy is implicit in this scenario. Following [23], 
for this decay scheme with the condition that the coherence is lost (baseline 
L » 1/Am 2 and the oscillatory part is absent), the flux of v a for flavour a 
on reaching the earth from a distant GRB at a distance L from the earth, is 
written as 

K(E) = J2J2€ b (E)\U bl \ 2 \U m \ 2 ex P (-AnL/(\ d ) l ) . (10) 

b i 

In the above, a, b are the flavour indices, i is the mass index and U is the MNS 
mixing matrix for neutrinos. The exponential term in the above equation is 
the decay term with (\d)i is the decay length given by, 

(A,), = ^kmJL^ (11) 

where E is the neutrino energy and «j(= mj/r), r being the rest frame decay 
lifetime and m 8 is the mass of the i th decaying neutrino. Thus, the decay 
lifetime in lab frame (~ Er/m) and hence the decay length has a strong 
dependence on the neutrino energy too. 
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The fluxes of u e , and v T from a distant GRB can now be calculated for 
this decay scenario using Eqs (1-4) and Eqs (10-11). Eq. (4) gives the total 
flux in absence of any oscillation or decay. Consdering the fact that at source 
the flux ratio of u e , and v T is 1:2:0, in absence of decay or oscillation, the 
u e flux (0*J and flux (</>l ) and v T flux (<f>l T ) can be written in terms of 

i _ dN° hB 
^ v ~ dE° bs ' 

c e = \<Pu , = \$ v = 2<p s Ue , <p: t = o (12) 

Applying the above equation for u e flux (f) Ue on arrival at earth and with the 
condition that only \ui) is stable, we get, 

^ = ^S E )\ U el?\Uel\ 2 + W u S E )WMUel\ 2 + 

<j)t e (E)\U e2 \ 2 \U e2 \ 2 exp(-47rL/(X d ) 2 ) + 
2<j ) s Ue (E)\U^\ 2 \U e2 \ 2 exp(-4nL/(\ d ) 2 ) + 
<Pt e (E)\U e3 \ 2 \U e3 \ 2 exp(-4nL/(X d ) 3 + 
2<p s Ue (E)\U^\ 2 \U e3 \ 2 exp(-47rL/(\ d ) 3 ) . (13) 

The first two terms on RHS of Eq. (13) (the terms without decay factor) can 
be written as 

<j>l e \U el \ 2 [\U el \ 2 + 2|fV| 2 ] = \U el \ 2 [l + \U^\ 2 - \U Tl \ 2 } (14) 

where the use has been made of the unitarity condition J2i UaiUu = S a b. For 
13 = (U e3 ~ 0), |C/ Mj | 2 - \U Tj \ 2 ~ 0. The other terms of Eq. (17) can also 
thus be simplified for 6\ 3 ~ 0. With similar approach for and v T fluxes (0„ M 
and <fi VT respectively) can be calculated. The expressions for three neutrino 
fluxes on arrival at the earth are simplified as (in terms of the source flux) 

0„ e = 0: e (E)[|f/ el | 2 + |f/ e2 | 2 exp(-47rL/(A ( i) 2 ) 
+ |£/ e3 | 2 exp(-47rL/(A d ) 3 ) 

0„ M = 0: e (E)[|f/ Ml | 2 + |^ 2 | 2 exp(-47rL/(A (J ) 2 ) 
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+ \U^\ 2 exp(-4nL/(X d ) 3 ) 



4> Vt = <p s Ue (E)i\U Tl \ 2 + \U T2 \ 2 exp(-4nL/(\ d ) 2 ) 

+ \U r3 \ 2 exp(-4nL/(X d ) 3 ) (15) 

From Eqs. (10,13,15) one can see that in case of L >> X d , the decay effect 
is washed out. The decay lifetime (r) becomes so small that the neutrino 
under consideration decays completely much before it reaches the earth. Un- 
der this condition the flavour ratio takes the form <p Ue : <p v ^ : <J) Ut = \U e i\ 2 : 
\U e 2\ 2 '■ \U e 3\ 2 (from Eq. 15) as discussed in [23]. In this case, the actual 
distance of UHE neutrino sources such as GRB is not very important and 
one can work with the diffuse GRB neutrino flux such as given in Ref. [22]. 
But, if Xd ~ L then one cannot neglect the exponential term in Eqs. (13) 
and (15) and knowledge of L is essential. In such decay event, it is useful to 
work with neutrino fluxes from GRBs with definite redshifts (z) (and hence 
different baselines L) in order to probe the effects of neutrino decay. In the 
present work, oscillation parameters are chosen as 6 23 = 45°, #13 = and 
012 = 32.31° and redshift z = 0.03 unless otherwise mentioned. 

It may be noted here, that for a GRB with redshift z = 0.03 (~ 10 21 km 
= L), the decay length A^ ~ 4 x 10 20 km for unparticle coupling A = 0.004, 
d — 1.3 and A w = 1 TeV at neutrino energy E u ~ 2 TeV. Therefore for 
neutrino decay into an unparticle, the neutrino decay length is indeed of the 
order of the baseline length L. 

2.2 Detection of UHE neutrinos 

The ultra high energy neutrinos from a GRB can be detected in a ter- 
restrial detector such as IceCube by detecting the secondary products like 
muons, tauons, electromagnetic or hadronic showers that are produced due 
to charged current (CC) and neutral current (NC) interactions of UHE neu- 
trinos with the terrestrial rock and detector material. The IceCube is a 1 km 3 
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ice/water Cerenkov detector in south pole ice. The secondary muons that 
are produced via CC interaction can be detected by track signal. The v T CC 
interactions produce the so called "Double Bang" events (track + shower) 
whereas v e CC interactions gives electromagnetic shower at the detector. In 
the present calculation however, we do not consider the "double bang" events 
as this is perhaps effective for a narrow energy range of 1 PeV to 20 PeV 
for kilometer cube detector [24]. For u T CC interactions we consider instead 
the decay channel (r — > v^v T ) of secondary tauons where the muons thus 
produced give track signals. The NC interactions of all neutrino flavours 
produce shower events. 

The total number of secondary muons induced by GRB neutrinos at a 
detector of unit area is given by (following [25, 26, 17]) 

r E° hs A ]\[obs 

S = jfj" dE° bs ^P SUIV (E?\ z )P»(E?\ £ thr ), (16) 

where P SU rv is the probability that a neutrino reaches the detector without 
being absorbed by the earth. This is a function of the neutrino-nucleon 
interaction length in the earth given by L int = {a tot (E° hs ) iV^} -1 , (N^ is the 
Avogadro number and a tot is the sum total of CC and NC interactions) and 
the effective path length X(8 Z ) (gm cm" 2 ) for incident neutrino zenith angle 
Qz- Psuw takes the form, 

Ps Ur v(£° bs A) = exp[-X(^)/L int ] = ex P [-X(9 z )a tot N A }. 

InEq. (16), 

P M (£° bs , E thl ) = N A a cc (R(E: bs ; E thl )) , (17) 

with Na is the Avogadro number and <r cc is the CC interaction. The 
average range of muon inside the rock is given by 

1 rl-E tYl JE v ^ CCV^obs \ 

(i?(£° bs ; E thr )) = / A/i?(£° bs (l - V), £thr) Y ,V) (18) 

Jo ay 
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where y = (E° hs — E^)/E° hs . The range R(E^, E t hi) for a muon of energy 
is given as 

JE thI {dE^/dX} (3 \a + (3E thl J 

The average lepton energy loss with energy per unit distance travelled is 
given by [25] 

@> = — <>*' (20) 

The values of a and (3 used in the present calculations are 

a = {2.033 + 0.077 ln[£ M (GeV)]} x 10~ 3 GeVcm 2 gm _1 

(3 = {2.033 + 0.077 ln[E M (GeV)] } x lO-Wgm- 1 (21) 

for E^ < 10 6 GeV [27] and 

a = 2.033 x lO^GeVcmVr 1 

(3 = 3.9 x 10" 6 cm 2 gm- 1 (22) 

otherwise [28]. 



3 Calculations and Results 

We define a ratio R of the track events and shower events at the kilometer 
scale detector considered here. 

Muon track events . , 

R = —77, 23 

Shower events 

The muon track events are therefore the total sum of the events from CC 
interactions and v T CC interactions (as described in the previous section). 
They are obtained using Eqs (16 - 22) with the flux in the expression for S 
(total number of secondary muons) of Eq. (16) suitably replaced by 0„ or 
4> Vt (Eq. 15) as the case may be. The neutrino CC and NC cross-sections 
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(a cc and cr NC ) are taken from Ref. [22]. The effective path length used in 
the expression for P surv (Eq. (16)) is written as X{9 Z ) = J p(r(6 z , £)d£ where 
p(r(8 z , £) is the matter density inside the earth at a distance r from the centre 
of the earth and £ is the neutrino path length with zenith angle 6 Z . The earth 
matter density is taken from Ref. [26] that follows from preliminary earth 
reference model (PREM). The calculations are made for GRB neutrino zenith 
angle 6 Z = 100.9°. 

The shower events are obtained using the equation 

r dN r 1 da j 

N sh = J dE^P^Ev) x j -—P hlt (E u ,y) , (24) 

where, a- 7 = a cc (for electromagnetic shower from v e charged current inter- 
actions) or <j NC as the case may be and P int is the probability that a shower 
produced by the neutrino interactions will be detected and is given by 

P sh = pN A a j L (25) 

where p is the density of the detector material and L is the length of the 
detector (L = 1 Km for IceCube). Here too, the flux dN v /dE u in Eq. (24) 
is to be replaced by 0„ M , 4> Vt or 4> Ue (Eq. 15) as the case may be. 

From Eqs. (8,9), we see that the unparticle decay of neutrinos depend on 
the scaling dimension du (non integral number) and the neutrino-unparticle 
coupling strength X l J. In what follows, both X l J and just A both will sig- 
nify the same coupling. In order study how a possible unparticle decay of 
neutrinos can affect the UHE neutrino signal from GRB, we vary these pa- 
rameters and the subsequent variation of the ratio R (Eq. 23) is calculated 
with the formalism discussed so far. The unparticle decay of neutrinos con- 
sidered here, \v 2 ) and \v^) are considered unstable and subject to undergo 
unparticle decay while only is stable. Therefore, in Eq. (8) we need 
the masses m 2 and m 3 for unparticle decay of respective neutrinos. In the 
present calculations the value of m 2 is estimated from m 2 = \J Am^, where 
Am| 2 = wij — m\ (normal hierarchy) and Am\ 2 = 2.5 x 10 -3 eV 2 (from 
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atmospheric neutrino oscillation) for the present work. The value of thus 
follows. The cutoff scale Au is taken to be 1 TeV. 

Fig. 1 shows the variation of R with coupling strength X l J for three 
different values of scaling dimension du for a neutrinos from a GRB with 
fixed redshift z = 0.03. This redshift corresponds to a baseline length of 
3.8 x 10 21 Km. The ratio R becomes insensitive to the variation of coupling 
for higher values of du- This phenomenon can be understood from Fig. 2 
where the variation of decay lifetime r (in terms of ^) with d u is shown. 
This plot clearly shows that ^ increases with the increase of du- From Fig. 
2, for example, ^ ~ 10 18 when d u = 1.3. For UHE neutrinos with an energy 
~ 10 6 GeV coming form a GRB at 100 Megaparsec distance (L ~ 10 21 Km), 
exponential decay term in Eq. (10) tends to 1 and thus the muon track to 
shower ratio (R) tends to that for mass flavour oscillation. The value of 
the ratio R in case of mass flavour oscillation scenario for the same set of 
oscillation parameters and baseline length of GRB neutrinos (z = 0.03) is 
calculated as -Rmass- flavour = 3.23 and for no oscillation, R is computed as 
-Rnoosc = 5.6 for the same GRB. 

Since the coupling X l J appears at the denominator of the expression for 
^ (Eq. 8), for higher values of coupling the decay effects may appear if for 
those values, the decay length A^ becomes ~ the baseline length L. This 
is also apparent in Fig. 1. One also observes from Fig. 1 that for suitable 
values of scale dimension du and the unparticle coupling srength A, the ratio 
R can differ from mass flavour value (without the decay) and also from no 
oscillation value to a considerable extent (in fact the value R can even be 
very close to 0). Therefore, for certain conditions of unparticle scenario, the 
unparticle decay of neutrinos may indeed be probed by observing the UHE 
signal at a kilometer scale detector like IceCube. 

In Fig. 3, the variations of the muon track to shower ratio R with du 
are shown for four different values (0.0001, 0.001, 0.01 and 0.1) of couplings 
X l J . For higher values of du and for lower values of A, decay effects vanish 
and the value of R tends to that for only mass-flavour oscillation. A fact 



14 



that can be understood from Eqs. 8-11 and also discussed above. But once 
again, Fig. 3 shows substantial varitaion of the results from those obtained 
considering only mass-flavour oscillation, for certain ranges of values of du 
and coupling \ l J . For example, for du = 1.2 and A = 0.001, the value of 

R ~ 1.27 a variation of more than 60% from the mass-flavour oscillation 

value. Again, to obtain similar effects for higher value of du, the value of A 
is also to be increased. But perhaps, the unparticle coupling with neutrino 
may not be very high and the scale dimension du is to remain within the 
value 1 < du < 2. This may be mentioned that in Ref. [12] the unparticle 
decay of neutrinos is used to study the constraints on du, A w etc. 

The muon yield and muon to shower ratio are also affected by the GRB 
neutrino flux. The GRB neutrino flux has a d~ 2 (z) dependence, (Eq. 4) 
where d(z) is the distance of a GRB with redshift z. The flux also ofcourse has 
a E~ 2 dependence (Eq. 1). In an event of the unparticle decay of neutrinos, 
the variation of muon to shower ratio (R) signal at a kilometer scale detector 
like IceCube for different GRBs at different redshift distances is investigated 
in Fig. 4. The variations of R with z are plotted for three different values 
of unparticle coupling (A = 0.0001,0.001,0.01). Similar variation for pure 
mass-flavour oscillation and for no oscillation or decay are also shown in Fig. 
4. From Fig. 4, one sees that to obtain a variation of ~ 50% — 60% for the 
value of R from that of the pure mass-flavour oscillation value, the coupling 
is to be large or redshift z (and hence GRB distance d(z)) also should be 
large enough. For the latter case, however, the neutrino flux from GRB will 
be reduced (Eq. 4) which may affect significant detection. 

4 Summary and Discussions 

We have considered here the recently proposed "unparticle stuff" associated 
with the possible existence of a nontrivial scale invariant sector with scale 
dimension du- Such "unparticle stuff" appears to be a nonintegral number du 



15 



of invisible particles. Unparticles may open up unexpected phenomenology 
and in the present work we consider the decay of neutrinos to unparticles. 
The consequences of such decay process is explored with the example of ultra 
high energy neutrinos from cosmological Gamma Ray Bursts. We calculate 
how the unparticle decay of neutrinos, if exists, affects the neutrino fluxes 
(of all active flavours) from GRB on arriving at earth. We then attempt to 
estimate how the effect a possible unparticle decay of neutrinos can change 
the detection yield of GRB neutrinos from only mass-flavour oscillation (or 
no oscillation) at a kilometer scale ice Cerenkov detector like IceCube. The 
GRB neutrinos are chosen because of the astronomical baselines of GRB 
neutrinos (for a terrestrial detector), the oscillatory part is averaged out 
and one obtains an overall suppression for muon neutrinos. For the present 
calculation we probe the ratio R of the total secondary muon tracks to the 
total shower events at the detector for various possible scenarios. The muons 
and the showers are obtained as the secondary particles following the CC 
and NC interactions of the cosmic neutrinos with the terrestrial rock and 
detector material. We obtain more than 60% — 65% deviation of the value of 
R from only the mass-flavour oscillation scenario for some cases of unparticle 
decay of neutrinos considered here. 

Such deviations can be probed by an efficient detector. The volume of 
IceCube is 1 Km 3 but this is the volume over which the optical modules 
(OM) are sparsed. Even if a CC interaction vertex occurs in the ice outside 
this volume, the resulting secondary particle can be detected by IceCube if 
it enters the specified detector volume. Thus the effective volume of the de- 
tector is in fact more than 1 Km 3 . Although, IceCube detector has better 
energy and angular resolutions, there are backgrounds from cosmic muons 
and atmospheric neutrinos. Simulation study of IceCube detector by Ahrens 
et al [29] indicates that the cosmic neutrino signal from the diffuse cosmic 
neutrino flux is well below the atmospheric limit. A diffuse cosmic isotropic 
flux results from the summation of the cosmic sources. But here, we are 
considering the flux from specific GRBs. Thus they are like point sources. In 
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this case, one expects an excess of events from a particular directions. This 
considerably reduces the background if the detector has very good angular 
resolution. The simulation done at Ref. [29] with 1° angular search cone 
shows an improvement of the measurement with the diffuse flux. Moreover, 
the atmospheric neutrino spectrum falls more steeply (~ E~ 3 ' 7 ) as against 
hard GRB neutrino spectrum (~ E~ 2 ) from the shock acceleration mecha- 
nism considered here. Considering all these, a very large deviation of the 
signal from the expected ones (e.g. the ones expected from mass-flavour os- 
cillation) may be probed with IceCube. But detailed simulation is required 
to estimate the deviation of the results that can be probed by the detector. 
In the present work, the purpose is to study the effects if neutrino decay 
to unparticles vis-a-vis undergoes mass-flavour oscillation and no oscillation 
and GRB neutrinos are considered for that purpose. 
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Fig. 2 Variation of r/m (eV~ 2 ) with scale dimension dn. Case shown for m = rri2- 
See text for details. 
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Fig. 4 Variation of ratio R with GRB distance (in terms of redshift z) for 
unparticle decay for three values of unparticle couplings with dy = 1.2. Also 
shown are the only mass-flavour oscillation case and no oscillation (no decay) 
case. 
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